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Abstract. Supernovae play an integral role in the feedback of processed material into the ISM of galaxies 
and are responsible for most of the chemical enrichment of the universe. The rate of supernovae can also 
reveal the star formation histories. Supernova rate is usually measured through the non-thermal radio 
continuum luminosity, but in this paper we establish a quantitative relationship between the [FeII]i 26 
luminosity and supernova rate in a sample of 1 1 near-by starburst galaxies. SINFONI data cubes are used 
to perform a pixel pixel analysis of this correlation. Using Bry equivalent width and luminosity as the only 
observational inputs into Starburst 99, the supernova rate is derived at each pixel and a map of supernova 
rate is created. This is then compared morphologically and quantitatively to [Fell] 1 26 luminosity map. We 
find a strong linear and morphological correlation between supernova rate and [Fell] 1 26 on a pixel-pixel 
basis: 

log ; = (1.01 ± 0.2) * log ; - 41.17 + 0.9 

yr l pc~ l ergs 1 pc 1 

The Starburst 99 derived supernova rates are also in good agreement with the radio derived supernova rates, 
which further demonstrates the strength of [Fell] as a tracer of supernova rate. With the strong correlation 
found in this sample of galaxies, we now qualitatively use [FeII]i 2 6 to derive supernova rate on either a 
pixel-pixel or integrated galactic basis. 



1. Introduction 

Supernova rate (SNrate) is typically estimated by the integrated non-thermal radio continuum emission. 
The connection between radio continuum emission and SNrate comes from the infrared-radio relation 
that shows a tight correlation between star formation rate and both radio and infrared emission JT], and 
references therein). In addition, near-infrared (NIR) observations of supernova remnants (SNR) often 
show strong [Fell] emission line flux coincident with the radio peak. Iron atoms are typically locked in 
dust grains and shock fronts of the SNR cause efficient grain destruction through thermal sputtering. This 
releases the iron into the gas-phase where it is singly ionized by the interstellar radiation field, making 
[Fell] a strong tracer of shocks, which are often caused by expanding SNR. Imaging and spectroscopy 
point to [Fell] emission being coincident with known radio SNR, indicating that [Fell] can provide at 
least a relative measure of SNrate. 



2. Sample 

We have selected 1 1 bright near-by starburst galaxies at distances ranging from 10-100 Mpc. A summary 
of the galaxies included in this study can be found in Table [TJ All galaxies are spiral galaxies ranging 
from SO-a to Sd with the exception of NGC 520, which is a merger galaxy. There are two edge-on 
galaxies in our sample, NGC 520 and NGC 3628. The far-infrared luminosities range from 1.4 x 10 10 L o 
to 3.2 x 10 12 L o , where Arp 220 is considered an Ultra Luminous Infrared Galaxy (ULIRG), with 
L > 10 12 L o , and NGC 6240, NGC 1614 and NGC 7552 are classified as Luminous Infrared Galaxies 



Table 1. Summary of galaxy parameters. 



Galaxy 


Morph. 


Activity 


Dist. 


i 


log L FIR 








J2000 




[Mpc] 


[°] 


[Lq] 






NGC 3628 


Sb pec edge-on 


HII LINER 


12.8 


847 ±2 


79.29 


10. 


14 


NGC 4536 


SAB(rs)bc 


HII Sbrst 


15.4 


1802±3 


58.9 


10 


,17 


NGC 1792 


SA(rs)bc 


HII 


13.2 


1210+5 


62.78 


10 


,22 


NGC 1084 


SA(s)c 


HII 


16.6 


1409±4 


46 


10 


,42 


NGC 1808 


(R)SAB(s)a 


HII 


12.3 


1001+4 


83.87 


10 


,55 


NGC 520 


pec 


Merger Sbrst 


30.5 


2162±4 


77.49 


10 


,81 


NGC 7552 


(R')SB(s)ab 


HII LINER 


22.5 


1611+6 


23.65 


11 


,03 


NGC 7632 


(R')SB(s)0° 


HII 


19.3 


1535+15 


82.44 


11 


,43 


NGC 1614 


SB(s)c pec 


Sbrst 


64.2 


4778±6 


41.79 


11 


,43 


NGC 6240 


SO-a 


LINER 


108.8 


7242+45 


73.0 


11 


,73 


Arp220 


Sd 


LINER Sy2 


82.9 


5420±6 


57 


12 


,50 



(LIRGs), with L > 10 L Q ,. Our sample excludes galaxies with very active nuclei, however we have four 
Low-Ionization Nuclear Emission-line Region (LINER) galaxies and one weak Seyfert 2 (Sy2) nucleus. 



3. Methods 

3.1. Extinction Correction 

In order to study the true emission line strengths, it is crucial to determine the amount of dust extinction 
obscuring each nucleus. The extinction is calculated at each pixel and an extinction map is created 
showing the regions most effected by dust. The Ay is determined using the extinction law at NIR 
wavelengths ofA A x A~ LS 0. 

To grasp the true morphology and magnitude of the line emission, extinction-corrected [Fell] i 26 and 
Bry linemaps are created using the Pa/?/Bry line ratio. The [Fell] map is used directly to find the [Fell] 
luminosity at each pixel. The Bry linemap, on the other hand, is used along with the K band continuum to 
calculate both the Bry equivalent width and the Bry luminosity, which are used as the only observational 
inputs to calculate SNrate from the Starburst 99 model, described in detail in Section l3T2l below. 

3.2. Calculating SNrate 

Starburst 99 (from here on referred to as SB99) is a tool that models spectrophotometric properties of 
star forming galaxies such as spectral energy distributions (SEDs), luminosities, equivalent widths, 
supernova rates and colors. It includes predictions of the variations in these properties as a function of 
starburst age. These models are calculated for two extreme star formation modes: continuous, where star 
formation proceeds continuously at a constant rate, and instantaneous (delta function) starburst. 

A near solar metallicity (Z = 0.02), Salpeter IMF (a = 2.35, M hw - lM Q ,M up = 100M Q ) 
and an instantaneous star formation mode are used. For the case of an instantaneous starburst, the 
burst is normalized to an initial starburst mass of 10 6 M©. Using the observed Bry equivalent width 
(EW(Bry)), which is independent of this normalization, the average age of the population dominating 
the emission in each pixel is calculated. From age, SB99 then provides the expected SNrate. Because 
of the normalization, the SNrate must be appropriately scaled by comparing the SB99 age-dependent 
prediction of ionizing photon flux (N(H )) to the observed ionizing photon flux, which scales with Bry 
luminosity. This comparison allows a scaling factor that is directly proportional to the initial mass and 



Table 2. Average equivalent width of Bry and the SB99 derived average ages and integrated SNrates 
using the instantaneous starburst model. 



Galaxy 


EW(Bry) 
A 


(Mvr) 


SNrate.w 


NGC 7632 


5.9 


7.0 


0.03 


NGC 3628 


7.5 


6.9 


0.01 


NGC 4536 


10.6 


6.7 


0.3 


NGC 1792 


4.1 


7.9 


0.007 


NGC 1084 


4.8 


7.3 


0.009 


NGC 1808 


7.2 


6.9 


0.06 


NGC 520 


11.7 


6.6 


0.3 


NGC 7552 


13.7 


6.7 


0.3 


NGC 7632 


5.9 


7.0 


0.03 


NGC 1614 


22.9 


6.4 


0.9 


NGC 6240 


3.6 


7.7 


3.6 


Arp220 


8.0 


6.8 


0.7 



initial star formation rate, allowing the SNrate to be scaled based on the true initial conditions of the 
region represented in each pixel. 

4. Analysis 

4. 1. Age and SNrate 

Table [2] presents the EW(Bry) followed by the SB99 results, including the average age and integrated 
SNrate using the instantaneous burst model for each galaxy. Although the analysis was done on a 
pixel-pixel basis, the numbers presented in Table [2] are averaged over the galaxy, in the case of age 
and equivalent width, and integrated over the galaxy in the case of SNrate. 

4.2. Qualitative Correlation 

To demonstrate the morphological correlation between the [Fell] and SNrate, a side-by-side comparison 
of the K band continuum, extinction corrected Bry flux, extinction corrected [Fell] i 26 flux, and SNrate 
is shown in Figure Q] for NGC 4536, NGC 1808, and NGC 7552. These galaxies were selected for 
presentation because of the high quality of the observations. Visually comparing the morphologies of the 
four different maps, we find that the SNrate most closely resembles the [Fell] emission. 

Focusing on NGC 4536, the Bry emission has a bright knot directly north of the nucleus with 
secondary emission peaks to the north-west and south-east of the nucleus. The [Fell] emission is 
concentrated around the galactic center with a long plume of emission coming from the nucleus and 
extending towards the north-west corner of the galaxy. There is a small knot of faint [Fell] emission at 
the location of the Bry peak. Similar to the [Fell] emission, the SNrate peaks around the center with an 
arm extending towards the north-east corner. 

NGC 1808 displays a clear starburst ring in Bry surrounding the bright nucleus. The [Fell] reveals a 
much less prominent ring configuration with more extended emission around the nucleus spanning the 
south-east to north-west plane. The SNrate map shows more ring structure than the [Fell] but it is quite 
diffuse and lacks the strength of the Bry peak in the north-east region of the ring. 

The face-on spiral, NGC 7552, also reveals a starburst ring, with little Bry emission in the nucleus. 
Focusing on the northern Bry knots, we see three 3 main knots, two slightly north-west of the nucleus, 
and one elongated knot directly north of the nucleus. Comparing the Bry emission to the [Fell] emission, 
we discuss only one of the two north-western Bry knots. In addition, the morphology in the elongated 
knot differs from that of [Fell] and appears as two distinct knots. As shown by the SNrate morphology, 




Figure 1. A comparison between the K-band continuum, Bry linemap, [Fell] 1.26 linemap, and the SNrate 
morphologies for NGC 4536, NGC 1808, and NGC 7552. 



the northern emission knots resemble the realtive flux and morphologies of the [Fell] emission much 
more closely than the Bry emission. Specifically, the top north-western knot is also missing in the 
SNrate map and the elongated northern Bry peak is resolved into two peaks, resembling the relative flux 
ratios seen in [Fell]. 

The comparison of the Bry and [Fell] emission to the morphology of the SNrate per pixel reveals a 
general similarity between the [Fell] emission and SNrate, a much closer correlation than the [Fell] and 
Bry or the SNrate and Bry. We emphasize again the only observational input into the SNrate calculation 
is the Bry luminosity and EW(Bry), the [Fell] flux is never used. However, the SNrate correlates more 
with the [Fell] emission than the Bry. This strongly supports [Fell] as a robust tracer of SNrate. 



4.3. Quantitative Correlation 

To test if the [Fell] to SNrate correlation holds quantitatively, we compare the [Fell] 1 26 luminosity to 
the SNrate, pixel by pixel. The linemaps were not additionally filtered with the exception of NGC 6240, 
which was filtered in regions where the Bry flux is insufficient to determine the true extinction. The pixel- 
pixel plot of SNrate against [Fell] luminosity for each galaxy is shown in Figure |2l Both the SNrate and 
[Fell] luminosity are scaled per square parsec such that each pixel represents the same physical area. A 
least squares regression was performed in log space combining the points from all the galaxies excluding 



NGC 1792, due to the poor signal-to-noise linemaps. This combined linear regression is plotted by a 
solid black line in Figure |2] for direct comparison and the relation is shown below. 

log VS " rat \ = (1.01 ± 0.2) * log [Fe " ]l 26 o - 41.17 ± 0.9 (1) 

yr v pc 1 ergs 1 pc 1 

The above equation is a power-law with a slope of nearly 1.0, which represents a linear relationship. 
The error is calculated from finding the slope and intercept for each individual galaxy, and finding 
the standard deviation of that spread. Therefore the error represents the variations between individual 
galactic fits. 

The main uncertainty in the SNrate calculation is caused by the observational errors in Bry, [Fell] 
and K band continuum. The calibration uncertainty dominates the error, this is 10% of the flux for all 
cases except NGC 1792, where the uncertainty is dominated by noise and is estimated to be 20% of the 
flux. If we alter the Bry flux by ±10%, this translates to a change in SNrate of 10%. In addition, a major 
systematic source of error is the choice of burst model. Since our calculation of SNrate depends heavily 
on SB99, we are limited to only two star formation models, instantaneous and continuous. In the next 
section, the dependency on the burst model is discussed in detail. 

5. Comparison to Continuous Star Fromation and Radio SNrate 

An initial assumption used to perform this analysis was that our galaxies are best represented by an 
instantaneous burst of star formation. The other extreme case, also modeled by SB99, is the continuous 
star formation scenario. It is unlikely that our galaxies are represented by either of these extreme cases, 
but rather somewhere in between a starburst that is extended in time. Since radio continuum is the 
classical tracer of SNrate, the radio flux densities provide an independent measurement of the accuracy 
of the SB99 derived SNrate. To estimate the radio SNrate, the equation given by 0], relating non-thermal 
radio luminosity and SNrate, is used along with VLA observations. The integration area is matched to 
that of SINFONI. Figure [3] compares the SB99 instantaneous SNrate to both the radio SNrate and the 
SB99 continuous SNrate. 

For all galaxies, except NGC 1792, NGC 4536, NGC 6240 and Arp 220, the radio SNrate closely 
matches those given by the SB99 instantaneous-burst model. As noted before, the NIR-emission line- 
maps of NGC 1792 are too noisy to be reliable. However, NGC 4536 is a nearly face-on galaxy with 
relatively low extinction and high signal-to-noise spectra, and the NIR data should be reliable. Thus, 
NGC 4536 is undergoing star formation that is closer to instantaneous but extended in time. Arp 220 and 
NGC 6240 have radio SNrates even higher than the SB99 instantaneous SNrates. Also, in Figure |2j the 
points for these galaxies lie above the best fit line, demonstrating a [Fell] excess. This excess is caused 
by merger-related shocks that are exciting the [Fell] in addition to the SNRs. These powerful merger 
(U)LIRGs may demonstrate a limit to the validity of the direct relationship between SNrate and [Fell] 
luminosity. In addition, determining an accurate extinction in (U)LIRGs presents a challenge. The under 
or over estimation of extinction leads to inaccurate [Fell] luminosities as well as Bry luminosities used to 
scale the SNrate. However, this can only be established by studying a larger sample of (U)LIRGs. In any 
case, it appears that the majority of (modest) starburst galaxies is well-represented by the assumption of 
a (nearly) instantaneous burst of star formation. In addition, the very good agreement between SNrates 
derived from the radio continuum, and from NIR data provides added confidence in the diagnostic 
strength of [Fell] as a tracer of SNrate. 

6. Conclusion 

Using SINFONI observations of 11 near-by galaxies, we performed a pixel-by-pixel analysis of the 
correlation between [Fell] i 26 emission and SNrate. First, Bry, Pa/3, [Fell] i 26, [Fell] i .64, and H2.2.12 
line fluxes were remeasured. The Bry line flux along with the K band continuum was used to find the 
equivalent width of Bry, which was then used as input into the SB99 model to estimate the age of each 




Figure 2. A pixel-pixel plot of SNrate, as derived 
from SB99, compared to [Fell] luminosity. Each 
galaxy is represented with a different color and 
the values are normalized to 1 square parsec. 
The black line represents the best fit powerlaw 
excluding NGC 1792. 



Figure 3. Comparison of the SB99 instantaneous 
SNrate (x-axis) to the radio SNrate (blue dia- 
mond), the SB99 continuous SNrate (red aster- 
isk) and for comparison the SB99 instantaneous 
SNrate (black asterisk). 



pixel. From age, SB99 provides a SNrate, normalized to an initial mass of 10 6 M o . The normalized 
SNrate is scaled by Bry luminosity to derive the true SNrate in each pixel. 

Comparison of the [Fell] i 26 luminosity to the SNrate (derived only from Bry equivalent width and 
luminosity) reveals a nearly linear correlation. This relationship is valid both on a pixel-pixel basis and 
for the integrated galaxy. For the integrated [Fell] luminosity and SNrate, the fit is remarkably tight with 
very little spread. However, to correctly determine the absolute SNrates, it is still critical to determine 
whether star formation has occurred in a (nearly) instantaneous burst, or has proceeded in a continuous 
fashion. SNrates derived from radio continuum observations may be used to distinguish these scenarios. 
Most of the modest starburst galaxies in our sample are best fitted assuming instantaneous star formation. 
However, we find that the relationship breaks down for (U)LIRGs in our sample. 

The use of [Fell] as a robust tracer of SNrate is a useful diagnostic tool. Specifically, it would allow 
measurements of SNrates from NIR observations of distant galaxies, where individual SNR can not 
be resolved. With the strong correlation found in this sample of galaxies, we can confidently say that 
[Fell] i 26 is quantitatively correlated to SNrate and can be used to derive SNrate on either a pixel -pixel 
or integrated galactic basis. 
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